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THREE- PARTICLE COLUSIONS IN A GAS OF HARD SPHERES 

W. R. Hoegy 

laboratory for Planetary Atmospheres 

and 

J. V. Sengers 

Institute for Molecular Physics 
ABSTRACT 

The triple collision integrals which determine the first density correction 
to the transport coefficients, are derived for a gas of hard spheres using the 
binary collision expansion. This expansion provides a convenient technique for 
classifying the contributions in terms of sequences of successive binary col- 
lisions between three molecules. Such sequences contain both interacting and 
non-interacting collisions. It is demonstrated that for three hard spheres all 
sequences terminate after four successive collisions independent of the inter- 
acting or non-interacting nature of the collisions. As a consequence the colli- 
sion integrals are related to a limited number of sequences with three and four 
collisions only. It is shown that equivalent results are obtained from the surface 
integral form of the triple collision operator, derived earlier by Green and 


Sengers. 
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THREE-PARTICLE COLLISIONS IN A GAS OF HARD SPHERES 


I. INTRODUCTION 

Non-equilibrium statistical mechanics predicts that the first density cor- 
rection to the transport coefficients of a gas can be represented by a term linear 
in the density n. That is, the thermal conductivity \ and the shear viscosity rj 
can be written as** 2 

k Tt k 0 + n + • * • 

V ~ V 0 + V x n + • ‘ * (1*1) 

As is well known, the transport coefficients \ 0 and v 0 in the low density limit 
are determined by the solution of the linearized Boltzmann equation. The 
Boltzmann equation takes into account only uncorrelated binary collisions. 

The transport coefficients k Q and 77 0 can be expressed in terms of collision 
integrals which involve the velocities of two molecules before and after a binary 
collision. 

In the past decade many investigators have derived the integral equation 
for the coefficients k x and 77 x of the first density correction. For a discussion ' 
and a bibliography we refer to a recent review article of Ernst, Haines and 
Dorfman.* This integral equation leads to collision integrals for the coefficients 
/Vj and v x that involve the effects of collisions between three molecules. For a 
gas of hard spheres it was shown by Green 3 and Sengers 4 that the triple colli- 
sion integrals can be reduced to a surface integral form, analogous to the binary 
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collision integrals determining k Q and t? 0 . That is, the integrals could be ex- 
pressed in terms of the initial and final velocities of three particles in specified 
sequences of successive correlated collisions. Thus the problem was completely 
reduced to a study of the dynamics of three spheres. 

In the earlier formulation, 3 ^ the dynamics of the molecules was described 
with the aid of a time-displacement or streaming operator S_ t (1 • • i). The 
operator S„ t (1 • • t) is a substitution operator which replaces the phase vari- 
ables of l particles at a given time by their values at a time t earlier. As an 
alternative approach several investigators have proposed to describe the dy- 
namics using a binary collision expansion. 3 

Here we reconsider the reduction of the triple collision integrals for a gas 
of hard spheres from the point of view of the binary collision expansion. This 
expansion enables us to decompose the triple collision integrals into a series 

'i 

of terms consisting of increasing numbers of successive binary collisions. The 
decomposition, thus obtained, differs in appearance from the decomposition 
derived earlier by the surface integral method.^ Nevertheless, we shall show 
in Section V that the two decompositions are completely equivalent. One of the 
advantages of the binary collision expansion is that certain properties of the 
triple collision operator, such as its symmetry, can be readily demonstrated. 
More importantly, we shall use the binary collision expansion formalism to 
eliminate those collision sequences that cannot occur according to the laws of 
mechanics. 
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The generalization of the Boltzmann equation to include the effect of three- 
particle collisions is obtained from a cluster ejqpansion of the Liouville operator. 
As a consequence the collision sequences involve not only interacting collisions, 
but also non-interacting collisions. In the latter type of collision the particles 
pass through one another, instead of exchanging their momenta. We shall prove 
that all sequences of binary collisions between three hard spheres terminate 
after four successive collisions. Several sequences of three and four successive 
collisions can also be ruled out as being physically impossible. Thus the triple 
collision integrals will be decomposed into a limited number of sequences with 
three and four collisions only. 

H. THE TRIPLE COLLISION OPERATOR 

A formal solution of the integral equation which describes the effect of triple 
collisions, yields for the coefficients and v^ 9 ^ 

K = *(pi) *' hAh) 

Vt = T5Et | d Pi®(Pi) : 1 3®(pi) • (2.1) 
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The functions A(p. ) and B(p i ) which are functions of the momentum p t only, 
represent the solutions of the linearized Boltzmann equations 


/p 2 5 \ Pi 

Mfo) = " V'5Sr _ T kT / nT 


I a B(p t ) = -*1) 


— ♦ o — ♦ 

Pi Pi 

m 


( 2 . 2 ) 


The mass of a molecule is denoted by m and <£(i ) is the normalized Maxwell 
distribution function 


<£(i) = (277 mkT)“ 3/2 exp(- p. 2 /2m kT) # 


(2.3) 


The tensor p® p t is the traceless tensor associated with the dyadic p t p 1# The 
operators I 2 and I 3 are linearized binary collision and triple collision integral 
operators. In the cluster expansion method developed by Green 8 and Cohen 9 
these operators can be related to the streaming operators S_ t 


f 2 

I 2 = I d* 2 0 12 A(12)*(l)tf>(2) £p n 

^ i* 1 


(2.4) 


I 3 = f dx 2 dx 3 0 12 {£(123 )-£( 12) 13) - **>(12) <*>(23) 


(2.5) 


3 

,£(12)} 0(1)0(2)£(3) 2^ Pii 


i ■ 1 



with 


l 

A(l--'fc) = limS. t (l ••■{) JJ S tt (i) . (2.6) 

i ■ 1 

We use the notation dx i to indicate an integration over the momentum p A and 
the position r , of particle i. The permutation operator P 1A interchanges the 
indices 1 and i. The operator 0* x j is a differential operator 




d * u u d 

dp A dfj <9Pj 


(2.7) 


where n. . is the pair potential of particles i and j . 

Equation (2.5) was the starting point of the previous reduction of the triple 
collision operator, 4 ' 6 However, for the binary collision expansion it is more 
convenient to use the € -method, first introduced by Zwanzig. 16 This method 
uses resolvent operators G(1 * »<i) which are the Laplace transforms of the 
streaming operators 


G(l- - 'fc) = l dt e~ €t S_ t (1 •• l) = [e+ «(1 •• 4)]“ 

Jo 


( 2 . 8 ) 


where U (1 • • l) is the Liouville operator for the l particles 


mi--*) = « 0 (i- 1) - 

l < i < j <£ 


(2.9) 
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and 


« 0<1 


l 

■L 


m 


dir 


( 2 . 10 ) 


i* I 


We also need the resolvent operator G 0 (1 • * l) f which generates the free stream- 
ing of the particles 


G 0 (1” *) = [e +«„<!•• *>]"». 

(2.11) 

The equivalence of the two methods to describe the first density correction 
was demonstrated by Ernst 9 Haines and Dor f man. 1 Using their notation 9 one 
finds 

I 2 - 1 im € B 2 ( € ) 

* e-0 

(2.12) 

I3 = lil J eB 3 ( € ) 
«-*0 

( 2 . 13 ) 

where 11 


€ Bj (e) = fdx 2 0 12 G(12)W(12)G; 1 (12)<£(1) #2) 

2 


i» 1 
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£ B 3 (e) 




{<9 12 G(12)(* la + * 23 ) G(123) W(123) 


+ e i2 G(12) g(12; 3) 

- & lv G(12)W(1.') (9 u 3(13) W(13) 

- 8 12 G(12)W(12)0 23 G(23) W(23)} 


3 

• G* 1 (123) <£(1)<6(2) <t>(3) Y' P u (2.15) 

£ sc 1 


The statistical factors W(1 • »£), 


W(l--^) = JI (l + f..) 


(2.16) 


are i*elated to the Mayer functions f . . 

f.j = exp(~ Uj j/'kT) - 1 . (2.17) 


The function g(12,3) is given by, 

g(12; 3) = W(i2)f 13 f 23 . (2.18) 

It is understood that the limit € ~» o is taken after all other operations have been 
performed. 
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The self-diffusion coefficient can be expressed in terms of the same oper- 
ator s, provided one deletes the permutation operators, P u . 

The first term of B 3 (e ) can be simplified using the relation 

£> 1j G(12)(^ 13 +0 2 ,)G(123) = ^ 12 G(123) - # 12 G(12) . (2.19) 

We find it convenient to symmetrize the operator B 3 (e ) 9 so that all three par- 
ticles play the same role. For this purpose we interchange the integration 
variables 2 and 3. Furthermore, we add the corresponding terms starting with 
0 23 , since they vanish upon integration oner p 2 and p 3 . Thus we obtain 

if 3 

«B 3 (e) = 7 I dx 2 dx 3 T(123, e)4>(l) <£(2)</>(3) P u (2.20) 

J i*I 

with 

T(123, e) = (^12 +< 9 i3 +< 9 2 s) G(123)W(123) - Y ^ a G ( a ) w ( a ) 

l I. t 

-HH e *i G (“>) { f *j * 9 *i ° 

The summations are to be taken over the three pairs 12, 13 and 23. 
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We express the 6 . . operators in terms of resolvent operators 


0 i2 + 0 i 3 + 0 23 = ^o” 1 (123) “ G” 1 (123) (2.22) 

& 12 = Gj* 1 (12) - G" 1 (12) = Gg 1 (123) - G _1 (12) (2.23) 


so that 


T( 123, e) 


G 0 _1 (G( 123) - G„) W(123) - G “ 1 (G(a)-G 0 ) W(a) 


“2121 G °' 1 ( G W ' G ») W ( a i){ f “ 2 +G 0 _1 ( G K)-°o) W K)} 


Go ” 1 


(2.24) 


From now on we oirat the arguments of the free streaming operator G 0 . 

It is our purpose to decompose this operator into a sum of operators each 
of which transforms the initial momenta of the particles to the momenta after 
a specific sequence of collisions. The equations quoted in this section apply to 
molecules interacting with any spherical symmetric short range repulsive pair 
potential. The remainder of this paper is specifically concerned with a gas of 
hard sphere molecules. 

in. BINARY COLLISION EXPANSION 

A detailed study of the binary collision expansion for the case of hard spheres 
was recently presented by Ernst, Dorfman, Hoegy and Van Leeuwen. 5 The 
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resolvent operators G(1 • •■{) can be related to binary collision operators by 


G( 1 •• -£)W(1 •• -t) 


g<,+g °Z] t “ g<i 


a 


W(1 •• i) 


(3.1) 


W(1 • • - 6 ) G(1 • • •£) = W(1 • > l) 


G 0 +G(l 


"*>L 


T G ft 

a o 


(3.2) 


While G(1 • • l) and W(1 • • l) do commute, the binary collision operators to be 
used depend on whether the overlap exclusion is written to the left or to the 
right* 

The binary collision operators T a and T x contain an interacting and a 
non-interacting term* 2 

T a = V + T a " (3.3) 

T a = T a i + T a n (3.4) 

To describe the effect of the binary collision operators we consider the param- 
eters that specify a collision between two hard spheres 1 and 2 with diameter a- 
(see Figure 1). An impact vector b x 2 is defined by 

^12 V 12 r i 2 * ^12 V 12 * ' ( 3 . 5 ) 

where "? 12 = fj - r 2 is the relative position, v 12 = - v 2 the relative velocity 

and v 12 the unit vector in the direction of v 12 . For b x 2 <a we define a perihelion 
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vector 


o*k “ b 12 + v 12 “ b ^ (3.6) 

and two contact times 

T ia = (r 12 • 9 ia ± V^-bi 2 2 )/v 12 • (3.7) 

The operators 2 and T 1 2 are defined by 5 

X 12 = 0 ( a " b 12) 8 ( T r2) R l2 ( 3 - 8 > 

TA = " ®(°‘" ,b ia) 8 ( r ia) < 3 - 9 ) 

V 2 = " e(°--.b la ) S(T+ a ) (3.10) 


where 0(x) = 1 for x > 0 and 0(x) = 0 for x < 0 and S(r) is the Dirac 8-function. 
The operator R 12 transforms the velocities v 1 and v 2 into the velocities v x ' and 
v 2 ' before the collision 13 


B l2 v 2 = v 2 ' = v 2 + v 12 kk . (3.11) 

A product of T operators and G 0 operators can be interpreted in terms of a 
collision sequence, when read from left to right. For example, Figure 2 shows 
the collision sequences associated with the four terms of T 12 G 0 T 13 . Since 
we consider backward streaming, the diagrams are to be read from top to bottom. 
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The operators T t 2 and are only different from zero, when the particles 1 
and 2 are in contact at the top of the diagram. The terms in 2 G 0 T 13 require 
that the conditions for two successive collisions be satisfied, such that t “ 2 <t ~ 3 . 
(Note that G 0 and G 0 Tj n 3 may include situations where 1 and 3 are 
colliding, while particles 1 and 2 are still overlapping.) For Tj n 2 and Tj n 3 the 
time ordering refers to the contact time r+ 2 and r t + 3 , respectively. The oper- 
ator product transforms the velocities of the particles at the top of the diagram 
' to the initial velocities at the bottom of the diagram. A minus sign is associated 
with each non-interacting collision. 

We mention some relationships between the T operators and the Mayer 
functions: 


f. T 0 = 0 , T a f a = 0 (3.12) 

and 

Ta'ffl = f /S T a - T a fg = f^Ta, (a^/3) (3.13) 

We shall also need the commutation relation 

Go f a - ^ G 0 = G 0 (T a -T a )G 0 = G 0 ( T a n - T a n ) G 0 (3.14) 

For a proof of these relations we refer to Ernst et. al. 5 

The binary collision expansion is generated by successive iteration of Equa- 
tions (3.1) and (3.2). For the two particle resolvent operator this procedure 
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terminates after one iteration 


G(a) W(a) = G 0 + G 0 T a G 0 (3.15) 

W(a)G(a) = G 0 + G 0 T a G 0 f (3.16) 


since*4 


T a G 0 T a = 0 , T fl G 0 T a = 0 , T a G 0 T a - 0 . (3.17) 

Substitution of (2.23) and (3.15) into the expression (2.14) for eB 2 (e) yields 

f 2 

e B 2 (e) = I dx 2 T, 2 0(l)^(2) ^P., , (3.18) 

i- 1 

which reduces to the familiar form of the Boltzmann collision operator. 

Similarly we substitute Equation (3.1) for l = 3 and Equation (3.15) into the 
expression (2.24) for the triple collision operator 

Z] TaG ° + 


r 


T(123, e) = 


L 


T G(123) W(123) 


L. a 




- 1 


(3.19) 


a l^ a 2 


As mentioned earlier, for a dynamical interpretation we read the terms from 
left to right. Therefore, we prefer to bring the overlap conditions W(123) and 
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f a to the left. Using the commutator (3.14) we obtain 


T( 123, e) 


7] \ K 123) G< 123) " 7] T - G <> 


+ 


l . a 


LL\.{ 


+ G 0 T a 


io ) 


Go' 1 


(3.20) 


V*2 


In order to express T(123, e) as a sum of terms that are convergent individually, 
we iterate (3.2) twice 


W( 123) G(123) 


W(123) 



T a G 0 + 


+ G(123)x 2:i: W* 


(3.21) 


a l? a 2 


so that 


T(123, e) 


2 ] 

a 



a l^ a 2 

(a 1 /^a 2 ) 


f /3T a G 0 T a 


+ 


+ 



T ai W(123)G(123)T a2 G 0 T a 


(3.22) 
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In Equation (3.22) we have indicated explicitly that the indices in two successive 
T operators refer to different pairs of particles. However, this condition is 
also satisfied automatically as a result of (3.17). Each term in (3.22) involves 
at least three conditions on the phases of the particles. In the first term, f^ f T a , 
two pairs of particles overlap; we refer to this term as the double overlap 
term. The double overlap term is the contribution according to the theory of 
Enskog: 15 it is the Boltzmann collision operator associated with one pair of 
particles multiplied with the excluded volume of the third particle.^ The second 
sum in (3.22) is a collection of terms that contain a single overlap condition. 
However, it should be remarked that the products of three T operators also in- 
clude single overlap configurations implicitly. 16 

The triple collision operator can be expanded into sequences with increasing 
numbers of successive correlated collisions by further iteration of (3.21). Noting 
that 

Ta, W(123)G 0 T a2 = (1+ f^)T ai G 0 T a2 < 3 ‘ 23 > 


we obtain 


K * £ ^ a i) 


T(123, £) 



T<*> (123, e) 


(3.24) 
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with 


T< 3 > (123, e) 


2 ] b f r T * 

a 

(a? fry? a) 



a x /a 2 


b T « 





+ 



a^aj/aj 

( a i^/3^ a 2 ) 


G ft T G 0 T a 

0 a j 0 


(3.25) 


T< 4 > (123, e) ■ ZZLL <** ^aj G 0 ^a 2 Gq Gq T a ^ ♦ (3.26) 

a 1“^ a 2^ a 3^ a 4 

( a i ^ fit a 2 ) 


T (s) (123, e) = UnZ ^gy.y.^y.T., 


27) 


a j/* a 2 / a 3 / ® 4 ^ a 5 
( a 1 ^ a 2 ) 


etc. 

In the derivation of (3.24) we have not used explicitly the fact that the operator 
T (123, €) operates on a function of the momenta alone. We shall show in Section 
V that with the latter restriction, T (3) (123, e) can also be written in the more 
compact form 


T< 3) (123, e) = 



a l^ a 2^ a 3 
( a a 2 ) 


( 1+ V) T * Go T a G 0 T e 


(3.28) 
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Each term T< s) (123, e) in the expansion (3.24) corresponds to sequences of s 
correlated binary collisions between the three particles. 


IV. REDUCTION OF TRIPLE COLLISION OPERATOR 

A decision as to when the expansion (3.24) terminates requires a study of 
the dynamics of three particles. We first mention some rules that are immediate 
consequences of the definition of the T operators. 


T aG 0 VG 0 T a - 0 

T aC 0 V G 0 T a = 0. 


(4.1) 


These equations express the fact that a pair of particles cannot recollide after 
a collision, unless the trajectory of at least one of the two particles is deflected 
by an interacting collision with the third particle. 



Go V, G 0 T a " 3 


G 0 T a 




°o T a " 2 G 0 T a " 3 


G 0 T a 



(4.2) 


The reason for (4.2) is that in the given sequence, none of the pairs is aimed 
to collide after the first three collisions. 

We list a number of lemmas, each of which expresses the impossibility of 
a specific collision sequence for three hard spheres. 


Lema X: 


: ai G 0 V 2 G 0 V 3 G 0 T o ‘ 4 G 0 T as = 0 


(4.3) 
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The lemma that three hard spheres of equal mass and diameter cannot undergo 

17 

more than four successive collisions was stated by Sandri and coworkers and 

18 19 

proved in detail by Murphy and Cohen, * However, in the original formulation of 
the lemma, no distinction was made between real collisions and interacting col- 
lisions. A real collision is a collision with the condition that it is not preceded 
by a non-interacting collision. 20 Strictly speaking, in order for three spheres 
not to undergo more than four collisions, it would be necessary to prove 

Ta t (l + Go V. ) G 0 T a 2 ( 1 + Go V, ) °0 V 3 ( 1 + G 0 V, ) G 0 V, 

x (l+G 0 T/ 4 )G 0 T as = 0 , (4.4) 

The factors 1 + G 0 ensure that the succeeding interacting collision is con- 
sidered only when it is real. Equation (4.3) implies (4,4), but not vice versa. 

The reason is that according to (4.4) the conditions for five successive inter- 
acting collisions conceivably could be met, in which case the sequence would 
be rendered hypothetical as a result of the interference of a non-interacting 
collision. An examination of Murphy’s proof shows that the latter argument is 
never used and that the conditions for five successive interacting collisions 
indeed cannot be satisfied regardless of whether the collisions are real or hy- 

<y\ 

pothetical. Thus Murphy’s proof justifies the stronger conclusion (4.3). 

In the formulation of the subsequent lemmas, we denote the three pairs of 
particles by a, j3 9 and 7. Thus we shall always use the convention a f (i f 7 r a. 
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Lemma 2: 


0 


(4.5) 


T a G 0 V S V G o T/3 3 
Lemma 3: T a G 0 Tg 1 G 0 T y l G 0 T a s 0 , (4.6) 

Lemmas 2 and 3 express the fact that those particular four collision sequences 
cannot occur.* 7 >^ Again we use the stronger interpretation for four succes- 
sive interacting collisions, independent of whether they are real or hypothetical, 
justified on the basis of Murphy 1 s proof. 


Lemma 4: 

T a G 0 Tg‘ G u T r " G 0 T a 

= 0 

(4.7) 


T a G 0 T/ G 0 f; G 0 T a 

* 0 , 

(4*8) 

Lemma 5: 

T„ G 0 T y n G 0 y G 0 T a 

= 0 

(4.9) 


T a G 0 V G 0 V G 0 T a 

“ o , 

(4.10) 

Lemma 6; 

f r T a G 0 T/G 0 T a = 

0 

(4.11) 


{ y \ G 0 V G 0 T a = 

0 , 

(4.12) 

Lemma 7: 

f y T a G 0 y G 0 T y = 

0 

(4.13) 


^7 G 0 Gq T y “ 

0 . 

(4.14) 


Lemmas 4-7 are new. They are various representations of a theorem which 
says that once the conditions for a recollision, T a G 0 T/ G 0 T a , are satisfied, 
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pair y cannot be in contact during the entire recollision process. A proof of 
the lemmas 4-7 is given in the Appendix. 

The lemmas presented above are to be supplemented with the following 
rules. First, the lemmas are valid when the left-most and right-most T operators 
are either interacting or non-interacting. Secondly, since a non-interacting col- 
lision does not change any of the velocities, the lemmas remain valid upon addi- 
tion of any number of non-interacting collisions. 

From these lemmas, we deduce the following theorem: 

T< B >(123, e) = 0 , for s>5. (4.15) 

Proof: 22 First we note that in each term of T< #) (123, e ), for s > 5, the subgroup 
of five left-most operators is equal to a term in T (5) (123, e). Thus it is suf- 
ficient to prove that each individual term of T (S) (123, e) vanishes. 

For this purpose we consider all possible combinations of five T operators. 
From lemma 1 we conclude that at least one of the intermediate collisions must 
be non-interacting. Equation (4.1) says that such a non-interacting collision 
cannot be inserted between two T operators with the same index. Equation (4.2) 
rules out the possibility that two successive intermediate collisions are both 
non-interacting. The remaining combinations are listed in Table I. Upon in- 
spection, we conclude that all terms vanish, since each contains a subgroup of 
four collisions that are ruled out by the lemmas. The appropriate lemmas for 
the individual terms are listed in the second column of Table I. The term 
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T a G 0 T^ n G 0 Ty G 0 Tg £ G 0 T y vanishes according to lemma 2, since T^ n can be re- 
placed with - G 0 Tg n . 

As a next step we investigate the terms of T (4) (123, e). The various com- 
binations of four T operators are listed in Table II. Again we conclude that 
several sequences are impossible. Thus (3.26) reduces to 


T (4) <123, O = ( 1 + f r)t T a G o V G o T « G » T r 

a//3 

+ • (4.16) 

Furthermore, the terms with t y can be deleted as a result of lemmas 6 and 7. 
In the terms of T< 3) (iaa. e) containing three T operators we can delete the 
overlap exclusion for the same reason. Summarizing our results we find 

T( 123, e) = T< 3 > (123, e) + T< 4) (123, e) , (4.17) 

with, 


T< 3 >(123, e) 


E 


(c&fryfia.) 


f T 

y a 


EE 

<x a 2 
(a 1 ?/3?a 2 ) 


f a T. G„ T, 


* m G 0 T ai + T ai G 0 G 0 Tg J , (4<18) 

a l? a 2 
C aj/yfl/a 2 ) 
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T (4 > (123, e) 



a l? a 2 
(a 


' ai GcV 2 G 0 T ai G 0 T /3 


G 0 T O! G 0 T/G 0 T l2 ] . (4.19) 

Thus the dynamics is restricted to a limited number of collision sequences with 
at most four successive collisions. 


V. COMPARISON WITH SURFACE INTEGRAL METHOD 

The fact that the triple collision operator can be decomposed into a sum 
of operators, each of which is related to a particular collision sequence, was 
demonstrated earlier by one of the authors.^ That derivation stated from the 
expression (2.5) for I 3 in terms of the streaming operators &(1 • • i). This 
operator was symmetrized following the same arguments as those used in the 
derivation of (2.20). For the configurational part of the integral we considered 
the positions of the three particles along their free trajectories. A surface in- 
tegral was obtained by integrating over r~ which is the time relative to the 
time of the first collision encountered when streaming backwards. The result 
was 


1 3 = 1 dX 2 dx 


J dx 2 dX 3 2Z H S [ T «\ ) Z] T „( a i : a 2 )<^(l)ci(2)^(3)2^P li .(5.1) 

*17*0,2 u i=l 
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where the summation over /x represents a summation over the six diagrams of 
Figure 3 # indicated by Rl t R2, Cl, C2, HI, and H2. Just as in Figure 2, the 
diagrams should be read from top to bottom. For convenience, we have re- 
tained the 8 -function to indicate that the integrand is evaluated at the time 
rj == 0 of the first collision. The successive collisions are time ordered ac- 
cording to the contact times r“ . Although not indicated explicitly, the trajectories 
in III and Cl should be continued until no further interacting collisions are en- 
countered. The effect of the operator associated with each diagram is to trans- 
form, the velocities at the top of the diagram into the velocities at the bottom of 
the diagram. In this convention each diagram of Fig. 3 represents actually two 
diagrams: one in which the third collision is interacting and another in which 
the third collision is non-interacting. The operator has a minus sign, when the 
number of non-interacting collisions is odd,, For a derivation of this surface integral 
form of the triple collision operator we refer to the earlier publications .4*6, 23 
It should be noted that the diagrams of Figure 3 do not specify the collision 
sequences completely, but that in addition some auxiliary conditions have to be 
imposed. The first auxiliary condition requires that the phases of the particles 

at the surface r~ = 0 should be restricted not only to non-overlapping con- 

i 

figurations, but also to receding phases. Thus any collisions that might occur 
when the trajectories are extended into the future, should be excluded. The 
second auxiliary condition says that all collisions up to and including the first 
non-interacting collision should be real . Lastly* as mentioned earlier, the three 
collisions of R1 and Cl could be followed by a fourth interacting collision. 
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with 


On comparing (5.1) with (2.20) we see that limT(123, e) should be identified 

e- 0 


K'23) = (5.2) 

a l/ a 2 M 

The six diagrams of Figure 3 represent a decomposition of the triple collision 
operator which differs in appearance from the decomposition (4.17) derived from 
the binary collision expansion. It can be shown that the auxiliary conditions men- 
tioned above represent a concise formulation of the combined effect of all operators 
T (s) (123, €) in (3.24) for s > 3, regardless of the validity of the lemmas quoted in 
the previous section. For convenience, we shall demonstrate the equivalence by 
considering only those collision sequences that are dynamically possible. Thus, 
we first investigate how the lemmas simplify the (a^ ; a 2 ) operators. 

To specify the auxiliary conditions explicitly, we introduce the following 
functions 


A* 

= e{cr-b a ) ) 

(5.3) 

Za 

= e(cr-b a )e(-T + ) 

(5.4) 

N a 

= e(b a -or) . 

(5.5) 

Thus A a = 1 when pair a is aimed to collide in the past, Z a = 1 when pair a is 
aimed to collide in the future, and N a = 1 when pair a is not aimed to collide in 
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either direction. These functions were used previously by one of us in a discus- 
sion of the Lorentz gas. 2 ^ Since f a can be written as 

f a = - e(o--b a )e(-r;)e(r + ) , (5.6) 

we note 

A a + Za + N a - f a = 1 . (5.7) 

The ( a l ; a 2 ) operators can be transcribed in terms of these functions 
and the binary collision operators used in the previous sections. As an example 
we consider T R1 (12; 13). 


S K 2 ) t ri< 12 ; 13 > = (A 1 3 + N 13 )(A 2 3+N 23 )T 1 i 2 (l+G 0 T 2 n 3)G 0 

x ^13 (* + ®0^23) (^12 + ^12 ^0^23) • 

The factor (A y +N y ) guarantees that the receding phase condition is satisfied 
for 7. The function 8 (r ~ 2 ) is incorporated in the first T 1 2 operator. The fac- 
tors 1 + G 0 T 2 n 3 guarantee that the second and third collisions are real. The last 
term T 12 G 0 T 23 gives the correction if a fourth collision is encountered. Lemma 
2 implies that such a collision can only involve p^. *es 2 and 3. Again it is 
understood that the limit e - 0 is taken just as in the preceding paragraphs. 
Similarly the other T (12; 13) operators ca*. represented by 

M T 1 2 ) ~ (^23 + N 23 ) Ti° 2 Gq T 13 Gq T 12 » (5.9) 
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( A i3 + N 13 )(A 23 + N 23 )T 1 i 2 (H-G 0 T 2 " 3 )G 0 



) T C1 

(12; 13) 

= ( A i3 + N 13 )(A 23 +N 23 )T 1 i 2 (l + G 0 T 2 " 3 )G 0 





X ^13 ( * + Go ( 3 + ^2*3 Go Tj 3 ) , 

(5.10) 

8 ( T fa) 

1 T C2 

(12; 13) 

( a 23 + N 23 ) G 0 T 13 G 0 T 23 , 

(5.11) 

S K 2 ) 

T hi 

(12; 13) 

= ( A l3 + N 13 )(A 23 +N 23 )T l i 2 G 0 T 1 " 3 G 0 T 23 , 

(5.12) 

s Ka) 

* ^H2 

(12; 13) 

a T n r. T n n T 

A 1 2 ^0 1 13 u 0 a 23 * 

(5.13) 


jin describing T R1 and T C1 we have used the fact that only one particular col- 
lision could be added to the three successive interacting collisions. Lemmas 4 
and 5 lead to further simplifications in T R1 and T cl 

M T 1 2 ) ^R1 ( “ (^13 + ^13 )(^23 + ^ 23 ) ^12 G 0 ^13 G 0 (^12 + ^12 G 0 ^23^®^^ 

0 ( T 1 2 ) ^C1 ~ (^1 3 + ^1 3 ) (^23 + ^23 ) ^*0 [^l l 3 G 0 ^23 

+ ^2 3 G 0 ^1 3 G 0 ^2 3 + ^13 G 0 ^1 2 G 0 ^2 3 + 3 G 0 ^2 3 G 0 ^ 13 ] (5,15) 

Note that the term T 12 G 0 T 2 n 3 G 0 G 0 T 2 l 3 G 0 T 13 in (5.10) vanishes according 
to lemmp. 2, since T 2 n 3 can be replaced with ~ T 2 n 3 G 0 T 2 n 3 , 

In contrast to the derivation of (4.17), the surface integral (5.2) was derived 
under the explicit assumption that the operator operates on a function of momenta, 
as is the case for the triple collision integrals (2.1). Therefore the two decomposi- 
tions of the triple collision operator (4.17) and (5.2) will only yield identical 
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results in the spatially homogeneous case. Since 

dx 3 G 0 l r ( r 2 i ) r 3 i* P i ’ P 2 ’ P 3 ) ~ 0 (5.16) 

if F is a function of the relative positions r 2 1 and r 3 , we can use a simplified 
version of the commutator (3.14) 

(T a n -T a ")G 0 = f B . (5.17) 

Furthermore we note that 

(A 2 3 +N 23) T 12 = (1 + f 2 3 " Z 2 3 ) T 12 = (* + f 2 3 + Vs G o) *1 2 ’ (5.18) 

The replacement of - Z 23 T 12 by T 2 n 3 G 0 T 12 represents a shift of the surface 
from r ~ 2 = 0 to r 2 + 3 = o, which is again justified in the spatially homogeneous 
case. Similarly A 13 + N 13 can be replaced with 

^13 + ^13 ~ ^ ^^13 G 0 > (5.19) 

where the term with f 13 vanishes due to the presence of the succeeding T 13 
operator. By substituting (5.18) and (5.19) into the expressions for (12; 13) 
and working out tne products, we can express T(123) in terms of products of T 
operators. 25 However, many terms vanish, again as a result of the lemmas 
quoted in the preceding section. Since the arguments are precisely the same 


27 



as those used in the reduction of T< 4 > (123, e) and T< s) (123, e), we do not dis- 
cuss the intermediate steps, but simply state the result. 

T(123) = LL‘K)L T m (a,; a 2 ) = T< 3 > (123) + T< 4 > (123) ,(5.20) 

f a 2 M 

with 


T< 3 > (123) 



*{a 2 
Jfa 2 ) 


1 

m 


' G„ T 1 

a, u 0 ‘a 2 


G 0 T a. 


+ T 


a. 


G n T a 
0 a. 


G 0 


] 


(5.21) 


t< 4 >(123) = J2 22 T «, 


G n T 1 G. T„ 

0 0^2 0 a j 


G 0 T /J +T a 




G o V G 0 T 


V a a 

(a 1 */3?a 2 ) 


“J 


. (5.22) 


The terms of T< 3 > (123) are precisely the terms associated with the six diagrams 
of Figure 3, if the auxiliary conditions are disregarded. The terms of T< 4 > (123) 
which are identical to those in (4.19) represent the effects of the auxiliary 
conditions. 

In order to show that (123) is identical to (4.18) we heed to shift the 
surface from t~ = 0 to r* = 0, when the first collision is non-interacting;. 
Using (5.17) we obtain 


T< 3 >(!23) = 


L L [\ °« t - 

a l^ a 2 


2 G ° T “l 


+ T. 


a l G ° T “2 


Go T /3 


1 S — 1 


i’LL 

a jr 

(CL x *(&a 2 ) 


f T„ 


G 0 T^ 


(5.23) 
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Repeating the procedure once more for the last term, we find 


* 


= he 

a l^ a 2 

(a l ?/3? i a 2 ) (a|//3/a 2 ) 

In the double overlap term f f a we can replace with , since both 
operators reduce to the Boltzmann operator in the spatially homogeneous case. 
Therefore on comparing these results with (4.18) we confirm the identity of the 
two forms of the triple collision operator. 

The operator can be transformed into a form which elucidates more clearly 
its symmetry upon time reversal. For this purpose we rearrange the terms of 
T< 3 > (123) 





f_ T. 


GoT/J 


+ f. f_ T 




(5.24) 


t< 3 > (123) = r 


b { y T a + 





a \? a 2 


T G ft T 

a | 0 a 


* G„ T 

2 0 a i 


T a j { f a 2 Go + G 0 T a 2 Go}J . (5.25) 


where from (3.14) 


fa 2 G 0 + G 0 T a2 G 0 = G 0 f aj +O 0 T aa G 0 . 


(5.26) 
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Using (5,26) and lemma 7, T< 4 > (123) can be written 



The time reversed operator is obtained by reversing the order of the operators 
and interchanging T and T. As mentioned above, the double overlap term 
ip T a is itself symmetric when operating on a function of momenta. Thus 
T< 3) (123) is symmetric upon time reversal according to (5,25) and (5.26) and 
(123) according to (5.27). This time symmetry is not sufficient to prove 
a generalized H-theorem, but it does ensure that the matrix elements of the 
triple collision operator in a Sonine polynomial representation^ are symmetric. 


VI, SUMMARY 

Using the binary collision expansion we have decomposed the triple col- 
lision operator into a series of terms related to collision sequences involving 
increasing numbers of successive correlated binary collisions. The expansion 
turns out to be equivalent with a decomposition of the triple collision operator 
derived previously by a surface integral method. We have presented some dy- 
namical lemmas which imply that all collision sequences between three hard 
spheres terminate after four successive collisions regardless of whether the 
collisions are interacting or non-interacting. This is a generalization of a lemma 
presented by previous authors which stated that three equal spheres cannot 
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undergo more than four interacting collisions. The dynamical lemmas lead 
also to some interesting simplifications in the sequences of three and four suc- 
cessive collisions. 
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' APPENDIX 

Our proof of the new lemmas 4-7 is based on an analysis of the recollision 
sequence, G 0 G 0 T * 2f (see Figure 4). The contact times of the three 
collisions are indicated by r “ , r A + , r 2 , r 3 ~, r 3 + . For convenience we follow 
the motion in the forward time direction. Lemmas 4-7 are consequences of the 
following theorem 

' r 32 < T ) > a for T 1 <T< r* . 'A.l) 

' > * . •" ’ 

V • 

In view of the symmetry of the recollieion sequence, it is sufficient to prove 

v 32 (r) > cr for T 2 <T <T+ . 

The theorem is most easily proved by examining the actual trajectories of 

the three particles in three dimensional space. For this purpose we consider a 

coordinate system with the center of the action sphere of 1 at the origin O for 

times Tf S t < t 2 * The coordinate system is oriented such that for times r >r 
sp 

sphere 2 is in the XZ plane, moving in the positive Z direction (see Figure 5)/ 
The relative separation of pair 21 at time rf is the vector from O to A 

QA = ? 2 i ( T i + ) = (^cos^j, 0, cr sin & % ) 0 < < § • (A.2) 

The center of sphere 2 is at point A at t “ t + . For times V > r* it moves 
along line AB in the positive Z direction. 
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At time r = r 2 , the center of sphere 3 lies at point C on the action sphere 

of 1, 

PC ~ r 32 ( 7 * 2 ) ~ (“ o- cos sin 0 2 , cr sin <fi } cr cos cos & 2 j , (A.3) 

where 0 is the angle that OC makes with the XZ plane, and & 2 is the angle be- 
tween the Z axis and the projection of OC onto the XZ plane, 

hi a collision between two hard spheres the velocity components along the 
line of centers are exchanged. Therefore for times r > r 2 the center of sphere 1 
moves along the extension of line CO, i.e., from O toward D and the center of 
sphere 3 moves in the plane perpendicular to CO at C. 

We examine the conditions on the location of point C (the center of sphere 3 
at time r = r 2 ) such that pairs 12 and 32 aim to collide at some time r >r 2 . 
Since at time r - r 2 the center of 1 is at O and the center of 2 is between A 
and B, 1 aims to collide with 2 only when the path of 1 is directed toward the 
tangent plane perpendicular to OA at point A; this requires C to lie in the op- 
posite hemisphere from A. Now, since at time r = r 2 the center of 3 is at 
point C and the center of 2 is between A and B, and A and C lie in opposite 
hemispheres, it follows that r 32 (r 2 ) > cr. Therefore 3 can aim to collide with 
2 only when the plane of the path of 3 intersects line A3. The above conditions 
imply for 6 2 and <P 



Consider the distance from the points on the path of 2 to the plane of the 
path of 3. Since the plane of the path of 3 intersects AB 9 there are two points 
on line AB whose distance to the plane of the path of 3 is cr. We denote as point 
E the position of 2 at the earlier time r Ef when the distance |ef| = c r , For all 
times r<r E ,r 32 (r)>cr f so that the time r~ of first contact of pair 32 is 
greater than or equal to r E : 

r 4 ~ T E • < A - 5 > 

Next we consider the distance from the points on the path of 2 to the plane 
through the Y axis and line OD. Since the distance from point A to the plane is 
cr sin (d 2 ~ d l) <a , at some time r G > r x + the center of 2 is at point G t such that 
the distance from G to the plane of OD and the Y axis is equal to cr . Thus the 
time r* of the last contact of pair 12 must be smaller than or equal to r G , 

r 3 < r G ’ ( A - 6 > 

We examine the distance from point E to the plane of the Y axis and line 

OD. This distance is given by the line EO. Line EO lies in the XZ plane and is 

therefore perpendicular to the Y axis; it is also parallel to CF and hence per- 
pendicular to line OD. The length of EO is 

l EO l = cr cos ^^cos Q 2 » ( A *7) 


37 


and thus, according to (A.2) and (A. 4) EO >cr . This implies AE > AG and hence 


On comparing (A. 8) with (A. 5) and (A.6) we conclude 


An equivalent statement is 


(A.8) 


(A. 9) 


r 3 2 ( r ) > a . for 7*2 5 T 1 r 3 • (A. 10) 

Lemmas 4, 5, and 6, Equations (4.7) - (4.12), follow from (A.l) since they 
involve the recollision sequence a p a with pair y in contact or overlapping dur- 
ing the recollision sequence. Lemma 7 also follows from (A.l) since must 
be an interacting T operator in Equations (4.13) and (4.14), which results in the 
recollision sequence y fiy with contact of pair a during the sequence. 
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Table II 


Terms In T (4 > (123, e) 


Collision Sequence Lemma 

T. G 0 V G 0 T.‘ G 0 T, = 0 2 

T. G 0 V G 0 T.‘ G # T y / 0 
T a G 0 T,‘ G 0 T y * G 0 T a = 0 3 

f. G 0 V G 0 T y ‘ G # T 4 t 0 
T. G 0 V G 0 T." G 0 T y / 0 
T a G 0 T,‘ G 0 T y G 0 T. = 0 4 

T.G 0 T /) "G 0 T y ‘G 0 T a = 0 5 

f. G 0 T/ G 0 T y ‘ G 0 T, ? 0 
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Figure 3. The six collision sequences associated with the surface integral form of the triple col- 
lision operator .4/6/23 jj nes represent particle trajectories and the circles indicate collisions 

between the two particles whose trajectories are enclosed. The dotted circles Indicate non- 
interacting collisions. 
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